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We reformulate the dynamics of homogeneous cosmologies with a scalar field matter source with an arbitrary self- 
interaction potential in the language of jet bundles and extensions of vector fields. In this framework, the Bianchi — 
scalar field equations become subsets of the second Bianchi jet bundle, , and every Bianchi cosmology is naturally 
extended to live on a variety of . We are interested in the existence and behaviour of extensions of arbitrary 
Bianchi-Lie and variational vector fields acting on the Bianchi variety and accordingly we classify all such vector 
fields corresponding to both Bianchi classes A and B . We give examples of functions defined on Bianchi jet bundles 
which are constant along some Bianchi models (first integrals) and use these to find particular solutions in the 
Bianchi total space. We discuss how our approach could be used to shed new light to questions like isotropization 
and the nature of singularities of homogeneous cosmologies by examining the behaviour of the variational vector 
fields and also give rise to interesting questions about the 'evolution' and nature of the cosmological symmetries 
themselves. 



1. Introduction 

An Equivalence Problem asks whether two geometric 
objects (e.g., manifolds, metrics, differential equations, 
Lagrangians, cosmological models etc) are the same 
under a suitable change of variables. Symmetries of a 
geometric object are defined as self- equivalences of the 
object and the determination of the symmetry group 
of an object is a special case of the general equiva- 
lence problem. Two equivalent geometric objects have 
isomorphic symmetry groups and, indeed, symmetry 
plays a central role in equivalence since if the symme- 
try groups of two objects are not isomorphic (e.g., of 
different dimensionality), the objects cannot be equiv- 
alent. 

There are two main approaches to Symmetry The- 
ory or Equivalence Problems, that of Sophus Lie and 
that of EUie Cartan. Lie's approach originally at- 
tempted to classify all possible Lie groups of trans- 
formations on one- or two-dimensional manifolds and 
has recently attracted a lot of attention [Q, especially 
the newly significant role it plays in problems of dif- 
ferential equations and variational calculus. On the 
other hand, the coframe equivalence problem, which 
includes all other equivalence problems as particular 
cases, introduced and algorithmically solved by Car- 
tan, although it initially triggered a lot of activity 
(mainly due to the effort of Cartan's students), subse- 
quently declined in applications due to its calculational 
complexity. However, very recently it received renewed 
attention and was revitalized primarily by Olver (see 
[H and references therein). Lie's approach is more 
"analytic" than Cartan's, which is more geometric in 
nature and is implemented by using the full machinery 
of differential forms. 

Equivalence problems in gravitation are abundant. 
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For example, the conformal equivalence |^ between 
different higher order theories of gravitation and gen- 
eral relativity with an additional field is achieved via 
a Legendre transformation which is a particular case 
of contact symmetry. Usually the first step to solving 
an equivalence problem is the determination of sym- 
metries of the object in question. Since Einstein equa- 
tions are the Euler-Lagrange equations of the Hilbert 
action functional, many equivalent problems in rela- 
tivity are bound to involve variational (Noether) sym- 
metries, that is, point or generalized transformations 
which leave the action functional invariant. Since ev- 
ery variational symmetry of a (variational) problem 
is a (Lie) symmetry of the associated Euler-Lagrange 
equations, one may first determine the complete sym- 
metry group of the equations and then decide which of 
these are variational. 

A different kind of equivalence problem appearing 
in relativity is that of classifying the symmetries of par- 
ticular relativistic models. Although there is to date no 
systematic attempt to tackle such equivalence prob- 
lems, some interesting work has been done to deter- 
mine the symmetries of some of the Bianchi cosmolog- 
ical models. In particular, Capozzielo et al. classi- 
fied the variational symmetries of some of the Bianchi 
Lagrangians and gave (through the Noether theorem) 
some first integrals (in particular, those for the Bianchi 
types I and V). 

Since the Lie approach is more 'easily' implemented 
and is less geometric in nature, in this paper we focus 
on the implementation of this method to the problem 
of finding the Lie symmetry groups of all Bianchi cos- 
mologies in vacuum or with a scalar field possessing 
an arbitrary self-interaction potential. We also give, 
when possible, the associated first integrals, thus com- 
plementing the previous work Q. 

The plan of this paper is as follows. Sec. 2 gives 
a brief introduction to the basic ideas of the geomet- 
ric theory of symmetries which is fundamental to our 
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work. In Sec. 3 we apply the theory of Sec. 2 to ho- 
mogeneous cosmologies and, in particular, recast the 
Bianchi dynamics in the language of jet bundles and 
extensions of vector fields and give a full list of sym- 
metries for the Bianchi-scalar models with an arbitrary 
self- interaction potential. The use of these results is ex- 
emplified in Sec. 4 where we give explicit forms of first 
integrals and construct solutions in many new cases 
of Bianchi-scalar cosmologies. We conclude in Sec. 5 
by pointing out several possible directions for future 
research. 

2. Jet bundles, extensions and 
variational symmetries 

Although in this paper we focus on sets of ordinary 
differential equations, for the purpose of this section 
it is very convenient to consider a general set of par- 
tial differential equations involving p independent vari- 
ables X — (x^, . . . ,xP) and q dependent variables u — 
(m^, . . . , m') . The total space is the Euchdean spaceQ 
E = M X N = {{x, u) : X e M CRP, u£ N CR'^} C 
'^p+<! . We restrict attention only to the most elemen- 
tary type of symmetry that is, point symmetries which 
are defined to be local diffeomorphisms of E onto it- 
self, g : E — > E : (x, u) i-^ g{x, u) 

(x, u) = g{x, u) = {X{x, u), *(x, u)) (1) 

which act pointwise on E. The set of all such diffeo- 
morphisms forms a group, the symmetry group of the 
total space E . A basic example of a point transforma- 
tion is constructed by starting with a vector field on E 
(that is a section of TE) 

V = Cix,u)dxi + (f>"{x,u)du'', 

i = l,...,p, a=l,...,q (2) 

and considering its flow exp(tw) . This defines a 1- 
parameter group of point transformations of the total 
space E onto itself. (In what follows, in a slight abuse 
of language, we refer to a symmetry either as g or as 
its infinitesimal generator v.) 

A function / : M — > N, u = f{x) , («e a section of 
the bundle E) is invariant under a group of transfor- 
mations G if its graph F/ = {{x, f{x)) : x E M} C E 
is a G -invariant subset, which means that ^g E G 
and {x,u) € Tf, g{x,u) E Tf. is a regular p- 
dimensional submanifold of E although only trans- 
verse submanifolds of E give (locally) graphs of func- 
tions of the total space (see [Q, p. 107). We must dis- 
tinguish between an invariant function u = f{x) and 
an invariant of G: A function I{x,u) on E f\ such 

^In the language of vector bundles E is replaced by a vec- 
tor bundle over the base space M (see [Q) and M and N are 
differential manifolds, not necessarily submanifolds of Euclidean 
space. 

■^Notice that for this definition x and u are not necessarily 
dependent variables. 



that u(/) = 0. This is a direct generalisation, on func- 
tions defined on E' , of the usual concept of an invari- 
ant on M , f{x), which means a real- valued function 
/ e C°°{M) such that v{f) = where v = f (x)a^. . 
(In this case we have a first order, linear, homogeneous 
partial differential equation v[u\ = which is solved 
by the method of characteristics.) In all these cases G 
is the symmetry group and members of G are called 
the symmetries of the function. 

How do we characterize invariant functions? Since 
the graph Tf of an invariant function u = f{x) is 
determined by the vanishing of its components — 
f°'{x), a = l,...,q, a direct application of the in- 
finitesimal invariance criterion {cf. Q|, p. 65) gives, 
using (I), 

= v[u- - f-^:,)] ^ r{x,u) - e{x,u)^, 

a=l,...,q (3) 

for every infinitesimal generator v of the form (|^). 
We call the g -tuple in the right hand sides of (|^) the 
characteristic of the vector field (||) and denote it by 
Q"(x, /(^^). Therefore (^) states that a function u = 
f{x) defined on the total space E is invariant under a 
(connected) group of point transformations if and only 
if it satisfies the equations 

Q"(x,/(i)) = 0, a^l,...,q. (4) 

Since in what follows we shall be dealing with 
symmetries of differential equations and of variational 
problems, we need to have a definition of symmetry 
for functions that depend not only on the x 's and u 's 
but also on the derivatives of the dependent variables 

and the independent variables . This necessarily 
takes us away from the total space E with the coordi- 
nates (x, u) and into the higher-dimensional analogues 
of E called jet bundles. This will also require an ex- 
tension of the infinitesimal invariance criterion (|^) to 
deal with such functions. 

Therefore we need to develop some aspects (partic- 
ularly interesting from the point of view of the applica- 
tions in the next Sections) of the theory of jet bundles 
and extension of vector fields. Furthermore, we shall 
have to recast differential equations in a more geomet- 
ric language using the notion of a variety on the jet 
bundle. Below, we give a very rapid overview of only 
those elements of the relevant theory we need, referring 
to and for a more complete discussion including 
references. 

The notion of a jet bundle is very simple. We may 
think of it as a space whose coordinates are the p in- 
dependent variables x* , the q dependent variables 
and the derivatives of , a = 1, . . . , g, of orders 1 up 
to and including n. This is the n-th jet bundle of the 
total space E ^ M x N C W xR'^ which we denote by 
J" = J"^; = M X iV(") where M is the (base) space 
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of the independent variables and A^(") is the fiber con- 
taining the remaining variables (dependent plus their 
derivatives of orders 1 to n). (A more rigorous way is 
to define J" as the set of equivalence classes of C°° 
functions wherein two functions are equivalent at x if 
and only if they are in n-th order contact at x, mean- 
ing that their Taylor polynomials at x of order n are 
identical.) 

Having 'extended' the total space E to the jet bmi- 
dle J'^E, it is natural to define the n-th extension, 
/^"■^ , of a function f : M — > iV to be a section of 
J", ie , : M — > iV^") is defined by evaluating 

all partial derivatives of / of order 1 to n . The graph 
of the extended function , r^(„) = {(x, , 
will similarly be a p-dimensional submanifold of J" . 
Since the point transformation will also act on the 
derivatives of functions / : M — » N, we can de- 
fine the induced extended (point) transformation g^"^ ; 
J" — > J" by 



(5) 



on the n-th jet space. This will transform the graphs 



,(")r, 
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( 



of extended functions giving 'i 
p. 113). 

A smooth real-valued function F : J" — > R (or 
on an open subset of J" ) is called a differential func- 
tion of order n. An n-th order differential equation is 
defined by the vanishing of an n-th order differential 
function. The total derivative, D^iF , of a differential 
function of order n with respect to , is an (n + 1)- 
th order differential function defined in the usual way. 
For instance, in the case of one independent variable 
X and one dependent variable m, we have the follow- 
ing formula for the total derivative of F{x, u'^"^) with 
respect to x: 

F)xF — Fx ^ UxF^i -\- UxxEUx UxxxEUxx "t- • . . . (6) 

Obviously vector fields on E can also be extended to 
vector fields on J^'E . The n-th extension of a vector 
field V with the characteristic Q = (Q^, • . • , Q') has 
the following form (see for a proof): 

= C{x.u)dx^ + <^,/"(x,uW)5„,= , (7) 

where 0,/" = fjQ" + ^"7,^, with i = l,...,p, a = 
1, . . . , g, I J| = j = 0, . . . , n (J is an obvious multi- 
index) and summation over repeated indices is implied 
in the usual way. 

A point symmetry of a system of (partial) differen- 
tial equations is a point transformation g : E — > E 
with the property that, if u = f{x) , is a solution, then 
the transformed function u — f{x) is also a solution. 
Suppose now that we are given a family of differential 
functions Afs : J" — > R : (a;,u(")) A^(a;, «(")), in- 
dexed by P — I, . . . ,m. A set of differential equations 
of order n is defined by the simultaneous vanishing of 
a given family of differential functions: 

A^(x,J"))=0. (8) 



It is very important to view these equations as defining, 
or defined by, a variety 

Sa - {{x, u^"^) : Ap{x, ?/(")) = 0, /3 = 1, . . . , m}, (9) 

which is a subset of J" consisting of all points of J" 
which simultaneously satisfy Eqs. (^. Thus in our ge- 
ometric reformulation a set of differential equations is 
a (suh)set of some space. A solution of Eqs. (^) is a 
function u = f{x) on E such that the graph of its n*'' 
extension, r^(„) , lies entirely on Sa (this is translated 
traditionally as 'the function u = f{x) identically sat- 
isfies Eqs. (|)'. 

It follows that the defining property of a symmetry 
of a differential equation discussed above, as a trans- 
formation that maps solutions into solutions, can be 
geometrically reformulated in an elegant manner by 
simply requiring that 



g(")(^A) c Sa- 



(10) 



On the other hand, the infinitesimal invariance crite- 
rion, (3), immediately implies that the fundamental 
condition for v to be a symmetry of the fully regular^ 
system of differential equations ((|], p. 179) is 



„(") 



(A^) = 0. 



(11) 



Notice that w" acts only on those points on 5a which 
lie on a solution. Eq. ( ITTI ) typically results in a large, 
overdetermined, linear system of partial differential 
equations and its solution requires very lengthy calcu- 
lations which are usually performed using some of the 
computer algebra packages available (see below, also 

In (g) and cff depend only on the total space 
coordinates x and u, and this was assumed by Lie 
[||. Subsequently, he extended that dependence to in- 
clude first or higher order derivatives of the functions 
It" , a — 1, . . . , g. Such generalized vector fields pro- 
duce the so called contact and generalized symmetries. 
More recently, nonlocal symmetries 0, |^, containing 
integrals of the dependent variable(s), have found their 
use. Thus there is a hierarchical sequence of symme- 
tries: point, contact, generalized, nonlocal. 

Once the symmetries of a differential equation have 
been determined, one may effectively reduce the order 
of the equation by calculating first integrals. For a 
set of differential equations, A/3(a;,u("^) = 0, a first 
integral is a function /(x, w'-™-') which is constant along 
solutions, that is 



Dxl = 0. 



(12) 



This equation is typically a linear partial differential 
equation and the associated characteristic system pro- 
vides a means to calculate all first integrals. If all of 
them can be determined, (8) is said to be integrated. 



^This, roughly speaking, means that the system has solutions 
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The complete solution of the (nontrivial) charac- 
teristic system of (18) is rarely trivial and one seeks to 
ease the process of solution by imposing the condition 
that any first integral be associated with a symmetry 
of (8), i.e., we require that 

= 0, (13) 

in addition to (8). This equation is also a linear par- 
tial differential equation for I and so the number of 
characteristics is further reduced. Usually, however, 
the existence of suitable symmetries (not necessarily 
point ones) equal in number to the order of the sys- 
tem is required, but the cost is the reduction of one 
first integral for each requirement of invariance under 
symmetry imposed. 

The use of Lie symmetries to determine first in- 
tegrals can be extremely difficult since the calcula- 
tions to be performed after the symmetries are ob- 
tained are usually difficult. The fundamental theorem 
of Noether |^ provides an attractive alternative. A 
point (or generalized) transformation is called a vari- 
ational (or Noether) symmetry of the functional 

A[u]^ ( L{x,u^''^)dx, ncM (14) 

Jn 

if and only if the transformed functional agrees with 
the original one. The set of all variational symmetries 
forms a group and this variational symmetry group is 
a symmetry group of the associated Euler-Lagrange 
equations but not conversely (cf. d), p. 236). The in- 
finitesimal invariance criterion applied to a (connected) 
transformation group G gives the basic condition for 
G to be a variational symmetry group of A[u\ , namely, 

w(")(L) + Ldive = 0, (15) 

for every infinitesimal generator v , where L is the La- 
grangian of the variational problem and ^ — (C) is the 
vector of basic components in (2) {cf. p. 236). 

As in the case of the Lie symmetries of differential 
equations, the variational symmetries and consequent 
variational first integrals impose no requirements on 
the coefficient functions ^ and (f> apart from differen- 
tiability. Variational symmetries can be point, general- 
ized (including contact) or nonlocal. The latter are not 
of practical use for a local Lagrangian unless the non- 
local terms cancel in such a way that the first integral 
is local. 

Given the ease of computation of variational inte- 
grals once the variational symmetries are known, one 
may wonder at the interest in the more difficult calcu- 
lation of first integrals using the Lie symmetries of the 
differential equation. For instance, provided the Hes- 
sian of the Lagrangian with respect to the first deriva- 
tives is nonsingular, the Euler-Lagrange system is 2n- 
dimensional and gives a regular Hamiltonian system. 
Each variational symmetry reduces the dimension of 



the system by two. If there exist n independent vari- 
ational integrals, there are n integrals in involution 
and the system is integrable according to Liouville's 
Theorem. 

Consequently, variational symmetries are very at- 
tractive. However, the calculation of variational sym- 
metries is a closed procedure only in the case of vari- 
ational point symmetries. In the case of generalized 
symmetries some ansatz must be made and there is 
always a possibility that some symmetry and so its in- 
tegral will be missed. (There are cases for which vari- 
ational symmetries are known to be generalized and 
yet the integral corresponds to a Lie point symme- 
try. Perhaps the best-known of these are the Laplace- 
Runge-Lenz components which correspond to gener- 
alized variational symmetries linear in the derivatives 
whereas the corresponding Lie symmetry is point [p^.) 
The fear of incompleteness in the knowledge of varia- 
tional symmetries in complex systems is a sufficient 
incentive to supplement the variational method with 
the Lie method. 

In this paper, we restrict ourselves to point symme- 
tries in both the variational and Lie approaches. (The 
calculation of other types of symmetry is extremely dif- 
ficult even for scalar ordinary differential equations and 
often does not give more information than obtained 
from the point symmetries.) We have remarked that 
the variational integrals follow easily once the varia- 
tional symmetries are known. We use the Lie method 
to supplement the variational approach. Indeed, there 
are some instances in which we can combine the two 
methods in the case that a variational integral is invari- 
ant under a Lie symmetry, which may or may not be 
the variational symmetry of the integral. The knowl- 
edge of this integral, both Lie and variational, can help 
one in solving the associated Lagrange system for the 
other integrals. 

Our strategy is now clear and will be described in 
detail in the next sections. We start by defining the 
Bianchi total space Eb to be the space including time 
plus the dependent variables, extend Eb to the sec- 
ond Bianchi jet bundle and consider extensions of 
all basic functions defined on the Bianchi total space 
(such as the Ricci scalar, etc.) to J^. In this way 
the field equations become conditions on Jg defining 
a Bianchi variety, Sa,b ■ This is a subset of the second 
Bianchi jet bundle containing all Bianchi models and 
consequently the structure of Sa,b is of fundamental 
importance to our work. The point symmetries anal- 
ysed in the next section have the property that, when 
extended, they leave Sa.b invariant. (This property is 
also shared by any contact vector fields although in this 
case (which do not consider in this paper) the Bianchi 
Lagrangian forms will not be invariants of the Bianchi 
jet bundle and will have to be modified to include Car- 
tan extensions known as coframes {cf. |2|).) 
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3. Point symmetries of homogeneous 
cosmologies 

In Bianchi homogeneous but anisotropic models the 
spacetime metric sphts so that the spatial (time- 
dependent) part is given by 

gab{t) = exp(2A) exp(-2/3)ab 

where A plays the role of a time (volume) parameter 
and /3 is a 3 X 3 symmetric, traceless matrix which can 
be written in a diagonal form with two independent 
components by introducing the two anisotropy param- 
eters (3+, 13- : 

Pab = diag (^f3+ , - 1/3+ + , - 1/3+ - ^/3_ 

The general Lagrangian leading to the full Bianchi- 
scalar dynamics has the form (see e.g. pd|) 



„3A 



R* + 6A' 



(16) 



where R* is the Ricci scalar playing the role of a 
potential term. The Bianchi total space is Eb — 
{(i, A,/3+,/3_,</))} C M^. Then the first and second 
Bianchi jet bundles are Jg = Eb x {(A, /3+, /32, 0)} C 



and Jl = J],y. {(A,/3+,/3_,</>)}c 



pl3 



with m — A and n — 2. Explicitly 



The Euler-Lagrange equations for (16) can be con- 
sidered as Eqs. 
(H) becomes 

3 / 'n 1 



1 



Pi 



d 



_e-3>— ( 



„3A 



Pi- 



R*) 
3/3iA 

3/?2A 



I OR* 



3 dPi 
IdR* 



3 dl32 
30A + V 



0, 
= 

= 

= 0. 



(17) 



For Bianchi class A models the Ricci scalar R* as 
a function on a four-dimensional hypersurface of the 
Bianchi total space Eb has the explicit form 



R* 



1 
2 

-2/3 



-2A 



■( 



-%/3/32 



^NiN2N3il + NiN2N3)., 



and for class B 
R* = 2a^ e"2A 



N2Ns\ , 



(18) 



(19) 



with 

P = 



3a2 - N2N3 



(7V2iV3/3i + v/-3a2iV2iV3/32) , (20) 



where a, Ni, N2, N3 are the usual classification con- 
stants. For the subsequent symmetry analysis we set 



u = e 



V = e 



Pi 



w = e 



V3/32 



(21) 



In what follows we list the relevant symmetries of dif- 
ferent Bianchi models in four cases according to differ- 
ent matter couplings: 

Case 1: Vacuum. It is obvious that the potential 
and (p terms are missing in the u equation and 
also there is no (j) equation. Hence we have a 
set of three equations with three dependent vari- 
ables u , V and w . This is the easiest case of all 
the models and, once having examined it, one is 
able to make inferences about the forms of the 
symmetries in the subsequent cases. 

Case 2: Scalar field with zero potential. In this 
case we have a set of four equations with four de- 
pendent variables u,v,w and 4>. The potential 
term is missing in the u and cf) equations. 

Case 3: Constant potential. The equations have 
the same form as in Case 2 plus a constant po- 
tential term in the u equation. 

Case 4: Arbitrary potential. It is the most gen- 
eral case. We have a set of four equations with 
four dependent variables. Moreover, in the u 
equation there is a potential function and in the 
(p equation we have a derivative of this potential 
with respect to 0. 

As we shall see, extra symmetries arise in cases 2-4. 
This is because the first three equations in (23) have (j) 
as an ignorable coordinate in cases 2 and 3 and so this 
is a symmetry of the fourth equation in (23). In case 
1, is not an argument of the coefficient functions. 

We also note that the actual calculation of Lie sym- 
metries in this paper use, to a certain extent, the pack- 
age LIE (cf. [^) which has been around for twenty 
years. The equation and type of symmetry sought are 
fed into an input file. The program computes the de- 
termining equations and then attempts to solve them. 
Usually, however, this proves impossible and in this 
case the operator can intervene manually. 

Bianchi Type I 

This is the easiest case since the Ricci scalar is zero. 

Case 1: The Noether symmetries are 

dt , vdtj , wdw , V log wdtj — 3w log vd^ . 

The additional Lie symmetries are tdt, utdu + ^t^dt 
and U(9,, . 
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Case 2: The extra variable (p gives three additional 
symmetries which are also Noether symmetries: 



dt, vdy, wdyj, d^, vlogw dy - 3w\ogv dw, 
v4>dv - I log vd^, w(j)dw - \ log wd^. 

The additional Lie symmetries are 

tdt, utdu + ^t^dt and udu- 

Case 3: The extra term of the constant potential 
V{(j)) = C affects the results and we obtain the fol- 
lowing Noether symmetries: 

dt, vdy, wdw, d^, t;logwa^ - Swlogua^, 

vcpdy - I logt;^^, w(pdw - | logwd^, 

and the Lie symmetries 

udu, c^^*9t + e^^*M9„, e^^'dt - e^^'udu. 

Case 4: Since this is the most general case, the num- 
ber of symmetries is obviously reduced. The Noether 
symmetries are 

dt, vdy, wdyj, V log wdy — 3w log vdyj, 

and the additional Lie symmetry is u9„ . 

Bianchi Type II 

The Ricci scalar takes the form R* = - c('^'3i"2A)/2^ 

Case 1: The system has the Noether symmetries dt 
and wdw ■ The additional Lie ones are tdt + udu and 

—2tdt + vdy. 

Case 2: We obtain the additional Noether sym- 
metry d^ and the additional Lie symmetry w4>dw 
+ ^ log wd^. 

Case 3: In this case the model has the Noether sym- 
metries dt, wdw and 9^, and the additional Lie sym- 
metry w(j)dyj + i log wi90. 

Case 4: In this case we obtain dt and wdy, as Noether 
symmetries and udu + ^vdy is a Lie symmetry. 

Bianchi Type III 

The Ricci scalar is R* = ge'^^ e^(^^=-^i). The 
symmetry dt is always a Noether one and 9^ is also 
Noether for the cases where it exists. 

Case 1: The Lie point symmetries are 

dt, tdt+^vdy, udu + ^wdyj, udu-^vdy. 

Case 2: In this case we obtain the following Lie point 
symmetries: 

dt, tdt+4:vdy, udu + ^wdyj, udu-^vdy, 

d^, v(j)dy + w(j)dyj - I (3 log w -I- log w) d^. 



Case 3: The Lie point symmetries are 

dt, udu + ^wdui, dij,, udu - ^vdy. 

Case 4: Arbitrary potential. We obtain the same 
symmetries as in the previous case apart from the d^ 
symmetry, which is lost. 

Bianchi Type V 

The Ricci scalar has the form R* = 6e^^^. 

Case 1: We obtain the Noether symmetries dt , vdy 
and wdw The additional Lie symmetries are tdt + udu 
and v log wdy — 3w log vdy, . 

Case 2: We find the additional Lie point symmetries 

2v(j)dy — 3 log vdif, , 2w(f)dyj — log wd^ 
plus the Noether symmetry d^. 
Case 3: We obtain the Noether symmetries 

dt, vdy, wdyj, d^. 
The additional Lie symmetries are 

v log wdy — 3w log vdyj , v<f)dy — I log vd^ , 

w(l)dy, - ilogwa^. 

Case 4: The additional Lie point symmetries are 

vdy, wdyj, vlnwdy — iwlnvdyj. 
There is only one Noether symmetry, dt ■ 

Bianchi Type VI, class A 

The Ricci scalar has the form 



R* 



. 1 p-2A 
2 



e4/3i + e-2(/3i-V3/32) ^ 2e''i+^'^2 



The only Noether symmetries are dt and dfj, in the 
cases where it exists. 

Case 1: The Lie symmetries are 

dt, tdt-\vdy-\wdu,, udu + \vdy + \wdu,. 

Case 2: There are the same Lie point symmetries as 

in the previous case. 

Case 3: We find the following Lie point symmetries: 
dt, udu + \vdy + ^wdyj, d^. 

Case 4: The system has the Lie symmetry udu + 
\vdy + ^wdy, . 
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Bianchi Type VI, class B 
The Ricci scalar takes the form 

R* = 2 (3a^ + 1) e"^^ exp 



3a2 + 1 



This model has only the standard Noether symmetries. 
If we apply Program LIE to this system, a number of 
difficulties appear. We set 



n = t;2/(3a^ + l) 



(22) 



and further simplify the system by defining the two 

constants 



S= i(3a2 + l), 



C = l/(3a^). 



(23) 



Case 1: In this case we obtain the following Lie point 
symmetries: 

dt, tdt + udu, 

- ^udu + ndn, \udu + ojO^. 

Case 2: The Lie point symmetries are 

dt , tdt + udu , - ^udu + ndn , 

n(j)dn + Lo<l)du, - iB^Clogw - fB^ logn d^. 



Case 3: In this case we find the following results: 
dt, -\udu + ndn, \udu + u}du,d^, 



ncpdn + uj<i)d^ - fS^ciogwa^ - lognS^. 



Case 4: We obtain the symmetries 

dt, ^udu + ndn, ^uOu+ojO^. 
The only Noether symmetries are the usual ones. 

Bianchi Type VII 

In this model the Ricci scalar takes the form 



R 



Case 1: There are three Lie point symmetries: 
dt, tdt - \vdy - fwc?^, udu + \vdv + | 

Case 2: We have the additional symmetry d^. 
Case 3: The symmetries of the system are 
dt, udu + \vdy + \wdu,, d^. 



Case 4: We lose the d^ symmetry of the previous 
case. In all cases the only Noether symmetries are the 
standard ones. 



Bianchi Type VIII 
The Ricci scalar is 



R* = -he 



-2X 



e^^i + e- 



-20 



1 (^e^0- + e-^l^-y 



Case 1: The symmetries of the system are dt and 
tdt + udu ■ 

Case 2: We obtain the additional symmetry Q^. 

Case 3: The system has only the two symmetries dt 
and d^. 

Case 4: The only symmetry is dt- 



Bianchi Type IX 
The Ricci scalar is 



R* = -^e-^^ 



,V302 _ „-\/3/32^ 



In Cases 1 and 4 the system has the single symmetry 
dt- The other two cases have the additional symmetry 

We now pass on to giving some examples of first in- 
tegrals associated with the symmetries obtained, post- 
poning their discussion till Section V. 

4. Applications 

In this section we use the general theory developed 
in Sec. 1 for obtaining first integrals corresponding to 
symmetries for Bianchi Types I, III and V and inte- 
grate to find some solutions. It is of interest that the 
expressions for the integrals are almost the same for 
Cases 2 and 3. Wc also describe how the existence of 
an additional variational symmetry restricts the possi- 
ble forms the potential can take to simple exponentials. 

4.1. Examples of first integrals 
4.1.1. Bianchi Type I 

Case 1: The results are given collectively in Table 1 
where 



1 

16 



{\oguf 



1/1 x2 1 , 1 

-l--(logu) +glog^+^ 



Using /3 , we find 



t-to 



I 



du 



(25) 
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Table 1: First integrals for the Bianchi model type I in 
Case 1 where there is no matter. 



Symmetry 


Integral 


vdv 

V log wdv — 3w log vdw 


11 — u^v/v 

12 ~ u^w/w 

u 



Table 2: First integrals for the Bianchi model type I in 



Symmetry 


Integral 


vdv 




Ii — u^v/v 


wdw 




I2 = u^w/w 








vlogw dv 




il^ 1 1 9 R 

u 4 6 3 




-~3w \ogv dvj 




— aarcsmh 



where 



(logw)'^ (logu)3 



(26) 



We formally invert (25) to find u{t) . The variables 
V and w follow from the integrals Ii and I2 , respec- 
tively. Hence we have 



exp / -^dt, w = exp / -^dt. 



We note that in this case we find the solution using 
only three integrals since they are separable Noetherian 
integrals. 



Case 2: Setting 
2 



C 



3VK 

2/3/4 
K ' 



K 



lit 
If 



K- 



1/2 



^2 



-'3 



we obtain the integrals listed in Table 2. (Note that /s 
may have arccosh instead of arcsinh, as we find after 
the calculation of the Lie integrals.) Integrating, we 
obtain the following solutions: 

u = [/3sinh(^-/5)-C]^/^ 
h 



exp 



rdt 



exp 



1 



X log 



/3tanh[(t 



C2 

'h)/2] 



p + ^/WTcp 



~Ctanh[(t - /5)/2] - /3 - y/(3^ + 



exp 



%dt 



1 



X los 



fJ2 

/3tanh[(t- /5)/2] 



/3+ VFTC2 



-Ctanh[(i - h)/2\ -(3^ ^/WTC^ 



Case 3: The integrals have the same form as in the 
previous case, the only difference being in the expres- 
sion for /3 where an extra term due to the potential is 
present, 

1/2 



2 /. 



C 



1^ (K-6P) - — 



2K 



(27) 



4.1.2. Bianchi Type III 

The integrals in Cases 2 and 3 are essentially the same. 
The integration becomes more involved in Case 3. 

Case 1: With the notations 
a — — 2 logu — — logw + — log 



f3 = VSlogu + ^\og 
v3 



(28) 



7 = log", 

the symmetry vdy + wdw transforms to dp with the 
associated Lagrange system 

dt da df3 d'y da d(3 d7 

The characteristics of this system are 

p — t, u = a, w — J, 

X ^ a, y = P, z = 7. 
Using these characteristics, we obtain 

dp du dw dx 
1 X z I (a;2 + t/2) _ 4e" 

dy dz 



(29) 



(30) 



-Zxy 

and thus 
du dy 

X 

du 

X '- 



(a;2 + j/2 + A^xz + 8z2) + | e" 



?>xy 



(31) 



da; 



(a;2 + ^2-) _ 4e« 



Ii = 



dp 
T 



du 

X 



2a^l'^ 
da 
a 



-/(a), 



(32) 
(33) 
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where 



f{a) = J e-6«da + Aja-i e"da. (34) 

(Note that P{t) follows from the quadrature of (|3l|).) 

The associated Lagrange's system for the symmetry 
dt is 



dt du dv d-y dx dy 

and gives the characteristics 

v = f3, 

y = P, 

From the system 

du dv dw dx 
X y z I (a;^ + y^) — 4e'' 
dy dz 



dz 

T 



(35) 



u — a, 
X — a. 



w = 7, 
z = 7. 



-3yx 



(a;2 +y^ + Axz + 8z^) + § e'" 



we have 
du 

X 



dz 



' I (2:2 + + 4a;z + 8z2) + I e" 



dz 3z2 3 



3 + y2 



2 e" 

3 X 



o- (36) 



wherein we already have x{u) and j/(u) from ( p^ ) and 
(^l|), respectively. Eq. ( p6| ) is a Riccati equation. On 
using the generalized Riccati transformation 



z = /w'(m)/w(u), 
we obtain the second-order equation 

CiJ CiJ LO X LO 

2 2 e" 



(37) 



3x' 



V 



3 a; 



(38) 



and, setting / ~ x/3 to remove the w'^ terms, we 
eventually obtain 



= u" 



9 x'^ +y'^ 
8 a;2 



2— ) 

a; 



(39) 



which is a linear second-order equation in w as a func- 
tion of u . So we deduce 7 as a function of a , 



du dui 

X z 

Hence 

t-tn = 



hi ■ 



7 



da. 



da 



, = / -^dt, 



(40) 



(41) 



so that this case is formally reduced to a quadrature 
and the solution of a linear second-order differential 
equation. 



Table 3: First integrals for Bianchi type V model in 
the first case where there is no matter 



Symmetries 


Integrals 


dt 


11 = u^v/v 

12 — u^w/w 


Vdy 


h = ]:{\oguf— - /i(u) 
2 u 



Case 2: In this case the integrals have the same form. 
The only difference is that 7 depends on 0. From the 
symmetry dt we have the extra integral 



which can be solved to find 
he-^^dt. 



(42) 



(43) 



Case 3: We obtain exactly the same integrals as in the 
previous case. We note that the second-order deriva- 
tive of 7 depends on the constant potential. 

4.1.3. Bianchi Type V 

We consider here only the problem of finding the Lie 
first integrals for Bianchi V in the first three cases, 
delaying the last case to the next subsection. In Case 
1 the integrals have the same expressions as those met 
in the Bianchi I (with a different form of /i and hence 
different gis). 



Case 1: The integrals are listed in Table 3 where 
^3 ^ 3(logu)2 31ogu 3 



(logu) 



(logu)^ + -(logu)^ 



Using the integral i we find 



t-tn 



-logu 
6 



du 



1 

36 



where 
9i{u) 



2u^l3 ^ 2uV(u) 
(logu)^ (logu)^ 



(44) 



(45) 



From (^^ we can deduce u{t) . The variables v and 
w follow from the integrals Ii and I2 , respectively. 
Hence we have 



exp 



h 



dt. 



w 



exp 



-dt. 



(46) 
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Table 4: First integrals for Bianchi type V model in 
the case where there is matter but no potential 



Symmetries 


Integrals 


dt 


11 = u^v/v 

12 — u^w/w 




h = 

h = i(logu)^- - f2{u) 
2 u 



(This solution uses only three integrals since they are 
separable.) 

Case 2: The integrals are shown in Table 4 where 



h{u) 



(logu) 



3(logM)" 31ogu 3 
2 2 4 



-{ill + ll + 2ll) 

{\oguf + ^(logw)^ + ilogit+ ^ 
2 6 36 

Using the integral 1^ , we see that 

du 



v/.92(w)' 



u-^. (47) 



(48) 



where 



92{u) 



2u'h , 2^2/2 (u)^ 



(l0gu)3 (logu) 

On the other hand, Ii and /2 give 
h 



V = exp 



rdi, 



w = exp 



rdt. 



Furthermore, using /a , we obtain 



exp / -|-dt. 



(49) 



(50) 



(51) 



We conclude that this case is reduced to quadratures. 

Case 3: The integrals have exactly the same form as in 
the previous case with a new function /3(u) replacing 
/2 (m) . We have 



1 



(log,)3 + ^(i2pL + ^ + 3 



1 



-C{\oguf--{Ml + ll + 2ll 



(logw)^ 



(52) 



4.2. Counteracting the symmetry breaking 
potential 

In this subsection we focus on Case 4. Note that, for 
instance, the vector field tdt is a symmetry for Cases 
1 and 2 but not for Cases 3 and 4 and this is obviously 
due to the potential couplings. In order to "fix" this 
problem we consider an additional term which "kills" 
any potential terms present in the A and cj) equations 
and becomes a symmetry for the other equations of the 
system. 

Consider the vector field 



V ~ tdt + add, 



(53) 



where a is a constant. Note that this is a symmetry 
for ([l7|b) and (p7|c) since the presence of the additional 
term does not affect these two equations. If we apply 
the new symmetry to (p7|a) and (p7|d), we expect to 
find restrictions on the form of the potential. 

Applying the second extension of (^ ) to (p7|a), we 
find that the potential must satisfy the constraint 



V" + -V' = 
a 

and similarly from (|r^d) we deduce that 

V' + -V = 0, 
a 

giving immediately 

V = Ke-^'t'l''. 



(54) 

(55) 
(56) 



Recall that the potential terms are responsible for re- 
ducing the number of symmetries in Case 2. We now 
consider the general point symmetry of Case 2 for 
which there is no potential, 

w = (A + Bt)dt + {C + Dt)dx 

+ {E + F4> + Gp2)dp, +{H + !((>- G(3i)dp, 



-102 d, 



(57) 



and apply its second extension to the set of equa- 
tions where the potential is unrestricted. Eq. (p^) 
implies that D ^ F = I = 0, which means that 
from the ten initially possible symmetries we lose three, 
the ones which correspond to the three vanishing con- 
stants. Also (Od) gives 



V = K 
= Ja + 2B, 



(58) 
(59) 



so that one obtains six independent constants as ex- 
pected (in Case 4 of this Type we found five indepen- 
dent symmetries). ???? 

A conclusion from this argument is that when one 
applies the general form of the symmetry, one ends 
up with the same symmetries as in Case 4 plus an 
additional one but with a restricted potential. 
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Our method can indeed be used in more general sit- 
uations. More general Bianchi cosmologies have extra 
terms in the Euler-Lagrange equations ( [1^ ) and one is 
forced to introduce a symmetry with an extra term. 
On considering the expression of the Ricci scalar in 
both classes and following the same reasoning as for 
Bianchi Type I, we are led to the symmetry 



tdt + add. + dx 



(60) 



If we take the second extension of (pOj), we conclude 
that for both classes A and B the potential is an ex- 
ponential function of 4> and has exactly the same form 
as in Bianchi Type I. 

5. Conclusions 

The geometric reformulation of homogeneous cosmolo- 
gies and the subsequent applications discussed in the 
paper allow some more general conclusions to be drawn 
regarding the dynamics of these cosmologies. All mod- 
els have in common the trivial Noether symmetry dt 
since all systems are autonomous. Case 2 of most of 
the models has the common symmetry . The sym- 
metries that one usually obtains in Cases 2 and 3 are 
combinations of the symmetries in the other two cases. 
Only Bianchi Type I in the case where the potential is 
constant gives a symmetry which involves the constant 
potential. This cannot occur in other models. Since 
Types VIII and IX are known to be the most complex, 
one does expect this to be reflected in the calculation of 
their symmetries groups and this was indeed the case 
as we obtained only the symmetries dt and 9^ . 

There are several directions in which one could 
extend research in this field. Firstly, the symmetry 
group calculations performed here can be extended 
to include other matter fields. Indeed, an interesting 
project could be to consider the effects of incorporat- 
ing a perfect fluid, electromagnetic fields, etc., on the 
number and nature of the symmetries found here. A 
detailed analysis of the nature and number of contact 
symmetries possibly present in homogeneous cosmolo- 
gies would be also very welcome. 

The well-known dynamical aspects of Bianchi mod- 
els must in some sense be reflected in the symmetry 
groups discussed here. How do our symmetries behave 
as these models expand? One does expect that neither 
their number nor their nature remains intact asymp- 
totically towards or away from singularities. Besides, 
some models are known to isotropise. Do their sym- 
metries evolve towards precisely those of the associated 
FRW models? We conjecture that they do, but this is 
a problem for the future. 
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